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We analyze the prospects for discovering and unraveling the nature of doubly- 
charged bileptons at a linear collider of the next generation running in its e~e~ 
mode. We stress the importance of initial state radiation, beam spread and polar- 
ization, and compute the discovery bounds. The gauge nature of vector bileptons 
can be determined by studying hard photon emission. 
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1 Introduction 

While the LHC offers an entry into the the high energy regime of the standard 
model with a significant opportunity for discovering new phenomena, a lin- 
ear electron collider of the next generation |J will provide a complementary 
program of experiments with unique opportunities for both discoveries and 
precision measurements. A major asset to fulfill this purpose is the versatility 
of linear colliders, as they can be operated in the four e + e~ , e~e~ , e~7 and 77 
modes, with highly polarized electron and photon beams. Moreover, starting 
from a centre of mass energy of several hundred GeV, it will later be possible 
to upgrade these machines into the TeV range. 

The e~e~ running mode @ is a particularly interesting feature of a high energy 
linear collider. One of the many promising processes which can be studied 
in this mode is the resonant production of doubly-charged bileptons. These 
particles are predicted by many extensions of the standard model, such as 
grand unified theories || , theories with enlarged Higgs sectors || , technicolour 
theories , theories of compositeness and theories which generate neutrino 
majorana masses @. They may appear as scalars, vectors or Yang-Mills fields 
and can be classified in a model-independent way || much like leptoquarks ||. 

The prospects for discovering doubly-charged bileptons through their s-channel 
resonance in e~e~ collisions have been considered previously within the frame- 
work of specific models [ |TT]|JTT[ ]p~|. We analyze these reactions here in more de- 
tails taking into account properties of the beams and effects of initial state 
radiation. We also study in a model-independent way the range of bilepton 
couplings and masses which can be probed in the e~e~ linear collider mode 
and we present ways of disentangling the different types of bileptons from each 
other. 

For this we present in the following section the most general classification 
of bileptons together with their interaction lagrangians with fermions and 
gauge fields. In Section 3 we study the discovery potential of bileptons and 
we demonstrate that the radiative return to resonance due to initial state 
radiation may enable the discovery of bileptons without having to resort to a 
painstaking scan of centre of mass energies. In Section 4 we investigate how 
to distinguish between scalar and vector bileptons, and to which extent a 
determination of the gauge nature of vector bileptons is possible. We conclude 
in Section 6. 



1 We do not agree with the cross sections (3.30) and (1) stated in the first two 
references of Ref. |i~Q|| , respectively. 
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2 Bilepton Classification 



We define bileptons to be bosons which couple to two leptons and which carry 
two units of lepton number. Their interactions need not necessarily conserve 
lepton flavour, but otherwise we demand the symmetries of the standard model 
to be respected. The most general renormalisable lagrangian of this kind is of 
dimension four and involves seven bilepton fields Lf, Lf + , LjL, L^, L®, Li 
and L3 + . It is given by |§ 



£ = - 


X\ 3 (I'PLVj-tjPM) 


+ 


X\ j Lf + I<P R £j 


+ 




+ 




+ 


V^Af Ll vtP L v 3 




Xl j Lt (itPtUi+llPLUi) 




V2X l i Li + llPdi 


+ 


h.c. , 



where the subscripts 1-3 label the dimension of the SU(2)l representation to 
which the bileptons belong and the indices i,j = e,fi,r stand for the lepton 
flavours. The chirality projection operators are defined as Pr % l = (1 ± 7s)/2. 

In the following we shall concentrate on the doubly-charged members of the 
scalar singlet L^~ , the vector doublet and the scalar triplet L^~ . The 
coupling matrices of the two types of scalars to different flavours of leptons A u 
are hermitian, while the coupling matrices for vectors are arbitrary. If CPT 
symmetry is imposed, the latter are either hermitian or anti-hermitian. 

For the sake of definiteness we shall focus from now on the universal diagonal 
coupling matrix 

A« = A PlJ , (2) 



which would be a natural possibility for gauge bileptons. It is straightforward 
to implement a different choice, which we will do at the end of this section. 
Within the framework of the coupling matrix (0), though, the negative results 
of muonium-antimuonium conversion searches [12| imply approximately the 



same upper limit for the three doubly-charged bileptons on the ratio X/tub, 
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where is the bilepton mass: 
A £ 0.5 m B /TeV . 



(3) 



Up to now we have made no difference between ordinary vectors and Yang- 
Mills fields for the doublet L 2fJ- bileptons. This distinction arises, however, 
when the interactions with the gauge fields are considered. In particular, the 
vector bilepton coupling to photon is described by the following lagrangian: 

£ = -1/2 {D^L 2u - D v L 2 rf (D^L 2u - D„L 2fl ) (4) 
- i k eQ L\L 2v (dy,A v - d u A fJ ) , 

where the covariant derivative is given by 

Df, = <9 M - ieQA^ , (5) 



and Q is the charge of the bilepton, i.e., Q = 2 for our purposes. If the vector 
bileptons are Yang-Mills fields, the parameter k in the second line of Eq. (§) 
takes the value one at tree level. On the other hand, for non-gauge vectors one 
expects a minimal coupling to the photon which is obtained when k vanishes. 



3 Discovering Bileptons 



Muons can be pair-produced in e~e~ collisions via the s-channel exchange of a 
doubly- charged scalar or vector bilepton, as depicted in Fig. [T]. The differential 
cross sections of these reactions are given by 



do-(Lr) (i + Pi)(i + p 2 ) a^ i 

dt ~ 4 2tt (s-m|) 2 + m|r| ' 

[RR] 

d<?{L 2 -) l-PiP 2 A 4 t 2 + u 2 



dt 2 8vrs 2 (s-m|) 2 + m 2 B T 2 B ' 

da(Lr) (1-Pi)(l-P 2 ) 2A 4 1 



(6) 



(7) 



dt 4 Ti (s-m|) 2 + m|r| ' 

" \LL] 



3 



where Pi and P2 are the polarizations of the incoming electron beams, mg and 
T# are respectively the mass and width of the bilepton, and A is the universal 
diagonal coupling (0). With this assumption on the structure of the leptonic 
couplings, the leptonic widths of the bileptons are given by []10| , |13| P1 





G = 


3, 




G = 


1 • 


L 3 


G = 


6 . 



(9) 



Taking into account the current limits (§) on the possible values of the coupling 
A, we see that bileptons as light as 200 GeV are bound to have a decay width 
which is significantly less than 1 GeV. 




Fig. 1. Lowest order Feynman diagram inducing the process e~e~ — > fj,~(j,~. The 
exchanged doubly-charged bilepton L can one of the scalar s L x ,L^~ or the 
vector Lo~ . 



In the case of the doublet and triplet bileptons, the different members of the 
multiplet need not necessarily have the same mass. If the doubly-charged mem- 
bers of the multiplets turn out to be heavier than the singly-charged ones, the 
non- leptonic decay mode L — > L~W~ can possibly also contribute to the 
total width However, only heavier bileptons can realistically accommo- 
date a mass splitting exceeding the mass of the W boson. We therefore ignore 
this possibility here and assume purely leptonic decays of the bileptons. 



Upon integrating Eqs 



the total cross sections are 



a 



s 



A J 



127T (s — m 2 B ) 2 + 'in 



2 p2 
B L B 



{In 


s = 


6 [RR] 




S = 


1 [LR] 


Us 


S = 


24 [LL] 



(10) 



If the decay width of the bilepton is significantly larger than the beam spread, 
it is a good approximation to assume monochromatic electron beams. In this 



We do not agree with the width stated in the text of the second reference of 
Ref. |0|. 
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case, the centre of the resonance can be precisely met, and the replacement 
s = m? B into Eq. (10) yields the cross section 



S 16vr 



a 



G 3m| 



This cross section does not depend on the coupling to leptons, since any 
decrease/increase of the strength of the interactions also results into a de- 
crease/increase of the decay width. In the observed cross section both effects 
result in a cancellation. Obviously, this result stops making sense for infinites- 
imal values of the coupling. Indeed, when the decay width is smaller than the 



beam spread, the latter needs to be taken in account ]TT]. If we approximate 
it by a boxfj]] of width R centrered around the resonance, the observed cross 
section results from the convolution 



m B +R/2 



a 



J (r B «R« rn B 



m B -R/2 



(12) 



where a typical value for the beam spread is [14 



R 



10" 



(13) 



Note that at a muon collider the corresponding spread would be of the order 
of 4 • 10~ 4 \/s From this fact alone, plus the luminosity enhancement at 



nominal energy associated with the much reduced Brems- and beamstrahlung 
rates of colliding muon beams, results the possibility of generating hundred 
times more doubly-charged bileptons at a facility. 

It may be that the resonance is so narrow that no significant tail extends on 
either side of the resonance. If this is the case and if the centre of mass energy 
were strictly confined within y/s ± R, it would be a long unpleasant task to 
find this resonance in the first place, through a minute energy scan. However, 
in practice initial state radiation and beamstrahlung will unavoidably perform 
a scan at lower energies. To obtain a rough idea of the effect, we approximate 
the centre of mass energy spectrum after initial state radiation to be of the 
form 



F(x,s) = — ln-4j 

7T 



1 



771g 1 — X 



(14) 



where a is the electromagnetic coupling strength, m e is the mass of the electron 
and x = s/s is the squared centre of mass energy fraction of the event. (The 



In reality it resembles more a saddle than a box [14]. 
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expression fll4|) can be derived, e.g., from the e~e" — > L 7 cross section 
(|20|), by singling out the large logarithm in the limit where <s.) Upon 
smearing the cross section ( |I0"D with the spectrum (0), we have 



1 

cr = y da; s) a{xs) 




5 si 

~ -2aA 2 ln — r (r fl < m B < v/i) . (15) 

G mi s — mi 



To estimate the discovery potential, we use the following scaling relation for 
the e~e~ luminosity 

C e - e - = 3.25- 10 7 s , (16) 



which closely corresponds to a luminosity of 25 fb _1 at y/s = 500 GeV which 
scales like the square of the centre of mass energy. This choice for the lumi- 
nosity is dictated by the latest e + e~ linear collider design reports |T3| and the 



fact that the e e mode will approximately suffer a 50% luminosity reduction 
because of the anti-pinch effect |1| . 



If we assume that observing one flavour violating e~e~ — > fi~ event already 
constitutes a discovery, we need an average number of — ln(l —p) Poisson dis- 
tributed events such that at least 1 event is observed with probability p. Hence, 
a predicted average of at least 3 events is needed to guarantee a discovery with 
95% confidence. 

These considerations imply that the minimal value of the lepton-bilepton cou- 
pling A needed to observe at least one e~e~ — > n~ event at the 95% confi- 
dence level on the resonance is according to Eqs ( |T2"||T3"|) 




A * 10 " 4 i- s ii * t ^ ■ <") 



Similarly, according to the cross section fllSD , far off the resonance the minimal 
coupling is 



A > 10~ 4 




\ 5 aln4 s 

\ mj 



Of course, the initial state radiation spectrum (0) is only an approximation, 
which at least has the virtue of yielding qualitative analytic results. We now 
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present the results of a more precise numerical convolution with a more real- 
istic resummed luminosity distribution [17|. We also assume a polarization of 
90% for the incoming electron beams. 



100 - 




500 



[GeV] 



Fig. 2. Energy distribution of the muon pairs in the reaction e~e~ — > [i~ fi~ 
at the centre of mass energy -y/s = 500 GeV. This processes is mediated by a 
doubly-charged vector bilepton of mass tub = 200 GeV which couples with 
strength A = 0.1 to the leptons (JJ). The peak at = tub = 200 GeV results from 
the radiative return to the bilepton resonance due to initial state radiation. 



In Fig. |2| we show the energy distribution of the muon pairs over its full 
range. The peak at 200 GeV corresponds to the mass of a vector bilepton L~^~ 
exchanged in the s-channel. The coupling to leptons is set to A = 0.1. Fig. ^ 
displays a zoom into this peak region for different values of the coupling. The 
effect of the radiative return to the resonance is important even for bilepton 
masses far below the nominal collider energy. This suggests that no painstaking 
energy scan may be needed to localize a possible narrow bilepton resonance. 

In Fig. H we show for several typical centre of mass energies the curves in 
the (me, A) parameter space which correspond to having a 95% probability of 
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Fig. 3. Same as Fig. || for different values of the lepton-bilepton coupling A and 
focusing on the muon energy region located around the peak. For the purpose of a 
clear display the two centremost bins for smaller couplings are depicted narrower 
than in the simulation. Their actual size is 1 GeV, like all other bins. 



seeing a e~e~ — > fi~fi~ event due to the s-channel exchange of a vector bilepton 
. As expected the sensitivity is maximal when mg = y/s and worsens 
dramatically below threshold. Beyond threshold, however, the sensitivity is at 
worst decreased by about an order magnitude. 

These results turn out to be actually quite close to the analytical predictions of 
Eqs (IrTyiSQ . As expected according to these relations, the corresponding limits 



for the singlet L x and triplet L 3 bileptons are shifted down by factors of 
a/2 and 2, respectively. 



While we have chosen the bilepton-lepton coupling matrix to be of the form 
(0), off-diagonal matrix elements (i ^ j) are in principle also allowed. For 
instance, in the case of scalar bileptons A u is a matrix of Yukawa couplings 
which, in general, need not be flavour diagonal and would lead to different 
final states in the processes e~e~ — > . In the presence of such non-diagonal 
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Fig. 4. Discovery prospects of a doubly-charged vector bilepton L^T as a function 
of its mass ms and coupling to leptons A (Q) for several values of the cm energy 
yfs. In the absence of any e - e~ — > event, the areas below and to the right of 

the curves cannot be excluded to better than 95% confidence. 



couplings the quantity plotted in Fig. | against bilepton mass m# becomes 



A ee A lJ 7y Sfc/(^ H ) 2 /3- It could thus be that the electron-bilepton coupling A ee 
is too small (much less than 10~ 4 ) to preclude the discovery of bileptons in 
e~e~ collisions. In this case they can still be pair- or singly-produced at a 
linear collider in e + e~, 77 and e~7 scattering WE . 



The low-energy constraints on the off-diagonal couplings ||[n§ are generally 
more stringent than the bound (|3]). Still, the combination of all linear collider 
operating modes will eventually probe the bileptons to exceedingly better 
sensitivities than any of the current experiments. 
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4 Determining the Nature of Bileptons 



A standard way of determining the spin of a particle is to examine the angular 
distribution of its decay products. Obviously, scalar decays must be isotropic 
whereas vector decays will display some non-trivial angular dependence. This 
is clearly verified by the differential cross sections 



However, since the availability of high electron beam polarizations (in excess of 
80%) is not a matter of debate JT^j, the polarization dependence of the scalar 
and vector total cross sections ( |T0"D will permit a much more straightforward 
determination of the spin of any discovered bilepton. Indeed, by flipping the 
helicity of one of the beams from the configuration with, say, two identical 
polarizations (Pi = P 2 = P) to the one with two opposite polarizations, (Pi = 
— P 2 = P) a scalar signal would be suppressed by a factor (1 — P 2 )/(l + P) 2 , 
whereas a vector signal would be enhanced by a factor (1 + P 2 )/(1 — P 2 ). Even 
with a polarization of only 80% (P = 0.8), this would mean a reduction or 
increase of the signal by respectively factors of almost 10 and 5. There is no 
way such a dramatic effect would escape detection. 

However, determining the gauge nature of a discovered doubly- charged vector 
bilepton is not such an easy task. It will necessarily involve a precise 
study of the interaction between photons and bileptons, according to Eq. 

We propose here a method based on the reaction 



e e 



L 



2 M 7 



(19) 



which proceeds through the three Feynman diagrams depicted in Fig. |5|. The 
third of these graphs involves the bilepton-photon coupling. As before, the 
lepton-flavour- violating decays of the produced bilepton guarantee an unmis- 
takable background-free signal. 



fwv7 



e~ — 



lyw N 7 




Fig. 5. Lowest order Feynman diagrams responsible for the reaction e e —* L 2 ^ 7 



The total cross section of the process (O) is given by 



10 



a 



aX 2 
3(s-m|) 



1 + P 2 
2 

1 _ P 2 



3A 2 



- (1 + A-A 2 ) 

+ 2k (1 + A - 2A 2 ) 

- 37 + 35A- 17A 2 

+ 6(A 2 -2A + 2) In 



(20) 



+ 



ml 



Although we have set the electron mass m e to zero wherever it is irrelevant 
in the final cross section (|20|), it is important to keep it finite throughout the 
calculation in order not to miss those terms which result from a cancellation 
of the electron mass in the numerator and the denominator. In particular, the 
helicity flip term on the first line is due to such a cancellation. As expected 
for such helicity flip terms, it vanishes at threshold where the photon energy 
is small. 



where 



m 



A = 1 — — and 



Close to threshold the cross section (^Cj) is dominated by the large logarithm in 
the penultimate term. Since this logarithm is solely due to the t- and w-channel 
electron exchanges, the gauge nature of the bilepton can only be revealed at 
higher energies. Away from threshold the last term, proportional to s/m 2 B , 
becomes increasingly more important. All the other terms are numerically 
irrelevant at all energies. 

At asymptotic energies the cross section ( p0|) saturates in the case of an ordi- 
nary vector with k — 0. In contrast, if the produced bilepton is a Yang-Mills 
field with k — 1, the cross section decreases like Although in principle 
this phenomenon may be used to distinguish these two cases, the onset of this 
asymptotic regime may be too slow to provide a useful test at the initial oper- 
ating energies of the linear collider. This is demonstrated in Fig. |6|, where we 
have plotted the energy dependence of the cross section (|20"D for an ordinary 
(k = 0) and a Yang-Mills (k = 1) bilepton of 200 GeV coupling with strength 
A = 0.1 to leptons. The difference is even smaller for heavier bileptons. 



11 




500 1000 1500 2000 



[GeV] 



Fig. 6. Cross section for the process e~e~ — > L^nl as a function of the centre of 
mass energy. The mass of the produced bilepton is 200 GeV and its coupling to 
leptons is 0.1. 



However, the gauge nature of the bileptons may also be determined from the 
angular distribution of the emitted photons. We do not give here the analytical 
form of the differential cross section corresponding to the process (|19l), since it 
is long and not particularly enlightening. However, we plot this distribution in 
Fig. ^ for an ordinary and a Yang- Mills vector bilepton of mass 200 GeV at a 
centre of mass energy of 500 GeV. The most salient feature is the existence of 
an radiation amplitude zero in the transverse direction for the gauge bilepton. 



To make efficient use of this difference in the angular distributions we per- 



form a Kolmogorov-Smirnov test | 20| . For this we first rewrite the angular 



distributions in term of the angular variable y = cos 2 9 in order to make ad- 
vantageously use of the forward-backward symmetry of this process. We then 
integrate these new distributions from y = up to an arbitrary value of y 
such that y < cos 2 8 m { n , where # min is the smallest acceptance of the detector. 
This defines two cumulative probabilities (for k = and k = 1) which increase 
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Fig. 7. Angular distribution of the process e e — > 7 for a bilepton of 200 GeV 
and a centre of mass energy of 500 GeV. 



monotonously as a function of y from at y min = up to 1 at y max = cos 2 9 min . 
Between these boundaries the cumulative probabilities for the two cases are 
different. Somewhere in the middle of the range the difference between the 
cumulative probabilities is maximum. The Kolmogorov-Smirnov statistic D is 
defined as the product of this number times the square root of the number of 
observed events. 

If we are dealing with, say, a Yang-Mills bilepton (k = 1), this cumulative 
probability is measured experimentally whereas the corresponding cumula- 
tive probability for the ordinary vector bilepton (k = 0) is computed. The 
Kolmogorov-Smirnov statistic is then given by 
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0.001 1 1 1 1 1 

500 1000 1500 2000 

m B [GeV] 

Fig. 8. Prospects for determining the gauge nature of a doubly-charged vector 
bilepton as a function of its mass tub and coupling to leptons A (^). The 

centre of mass energy increases from left to right as y/s = 0.5, 1, 1.5 and 2 TeV. 
The areas below and to the right of the curves do not allow a distinction between an 
ordinary vector and a Yang-Mills field to better than 99.9% confidence. The plain 
curves correspond to the case where the observed bilepton is an ordinary vector 
(k = 0) whereas the dashed curved correspond to the case where the observed 
bilepton is a Yang-Mills field (k = 1). 



where the luminosity £ is given by Eq. (|16|) and the indices and 1 refer 
to the cases k = and k = 1, respectively. If D = 1.95 the observed k = 1 
distribution has a probability of 99.9% to be different from a distribution which 
would have been observed if k = 0. Of course, if the observed bilepton is an 
ordinary vector, the corresponding Kolmogorov-Smirnov statistic is obtained 
from Eq. ([22] ) by interchanging the indices and 1. 

We have plotted in Fig. || the D = 1.95 boundaries in the (m^,A) plane 
for different values of the centre of mass energy. The sensitivity increases 
the further away the bileptons are tested from threshold. Far from threshold 
the sensitivity is slightly better if the observed bileptons are ordinary vectors 
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rather than gauge bileptons, because of the larger event rates. In this case the 
distinction between the two types of bileptons can be made for values of the 
coupling A as low as a few per mil. 



5 Conclusions 



We have investigated the potential of the e~e~ collision mode of a linear 
collider for discovering and studying doubly-charged bileptons in a model- 
independent way. For this we have used the most general renormalizable 
lepton-bilepton and photon-bilepton interactions for scalar and vector bilep- 
tons belonging to different representations of the SU(2) L gauge group, and 
we have studied the s-channel process e - e~ — > n~ n~ taking into account the 
effects of initial state radiation, beam spread and polarization. 

As expected, the linear collider sensitivity to bilepton properties is best when 
running on the resonance yfs = itlb- With the planned machine parameters 
couplings to fermions can be tested down to less than at the resonance. 
Far from the resonance there is still an important radiative return to the 
bilepton resonance due to initial state radiation, which makes it possible to 
discover bileptons for couplings down to less than 1CT 3 without having to 
resort to a collider energy scan. 

If a bilepton is discovered, its nature can be investigated with several methods. 
The simplest way to discriminate between scalar and vector bileptons is to 
make use of polarization, since flipping the helicity of the incoming electron 
beams would increase or decrease the production cross sections differently. The 
knowledge of its spin is thus a free by-product of the discovery of a bilepton. 

Recognizing Yang-Mills bileptons from minimally coupled vectors can be done 
by detecting hard photons in the process e~ e~ — > L^j. Because the angular 
distributions are different the gauge nature of a light bilepton can easily be 
tested for couplings as small as 1CT 2 or even smaller. 
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